Cold three-body collisions in hydrogen-hydrogen-alkali atomic system 
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We have studied hydrogen-hydrogen-alkali three-body systems in the adiabatic hypersplierical 
representation. For the spin-stretched case, there exists a single XH molecular state when X is one 
of the bosonic alkali atoms: ''Li, ^"^Na, ^^K, *^Rb and '^^^Cs. As a result, the only recombination 
process is the one that leads to formation of XH molecules, H-l-H-l-X— >-XHH-H, and such molecules 
will be stable against vibrational relaxation. We have calculated the collision rates for recombination 
and collision induced dissociation as well as the elastic cross-sections for H+XH collisions up to a 
temperature of 0.5 K, including the partial wave contributions from to 5~. We have also 

found that there is just one three-body bound state for such systems for and no bound 

states for higher angular momenta. 

PACS numbers: 



I. INTRODUCTION 

In the last decade, studies of three-body coUisional pro- 
cesses have attracted tremendous attention due to their 
great relevance for the rapidly growing field of cold and 
ultracold atomic gases [1]. In such systems, three-body 
recombination and elastic and inelastic atom-molecule 
collisions are of particular interest. 

Three-body recombination is a scattering process 
where three free particles collide, with two of them bind- 
ing to form a molecular state, converting the binding 
energy into the relative kinetic energy of the atom and 
molecule produced. Three-body recombination, there- 
fore, is important generically as it can shed light on 
binding in nuclear and chemical reactions. In ultracold 
atomic gas experiments, three-body recombination can 
lead to huge losses near a Feshbach resonance [2-4] and 
has been studied extensively to understand the lifetime 
and the stability of the gas samples [5, 7-11]. 

Elastic atom-molecule collisions are crucial for deter- 
mining the dynamics of ultracold atom-molecule mix- 
tures at the mean-field level, and inelastic atom-molecule 
collisions have a big impact on the lifetime of Fesh- 
bach molecules in such systems [12-14]. Furthermore, 
in the regime of large two-body s-wave scattering length 
a, achieved near a Feshbach resonance, three-body coUi- 
sional processes show universal scaling behavior with a 
as result of Efimov physics [6-9, 15, 16]. These univer- 
sal aspects have been observed experimentally in recent 
years [17-23] and confirm our understanding of three- 
body universal properties [16]. 

Ultracold three-body collisions, however, are only uni- 
versal when the dynamics arc predominantly determined 
by the long-range behavior of the atom-atom scattering 
wave function. When a system behaves universally, the 
complicated atom-atom interaction can thus be replaced 



with a much simpler model designed to reproduce the 
long-range wave function. It is in this context that the 
adiabatic hyperspherical representation has been applied 
to the calculation of three-body recombination. While 
it has proven very useful for getting deeper insights into 
this process, calculating recombination for chemically im- 
portant species using realistic interactions requires sub- 
stantial further technical development, but does not pose 
fundamental difficulties. 

Nevertheless, there are a few realistic systems that 
are sufficiently simple that recombination calculations 
are possible with the tools already available. For in- 
stance, recombination of helium atoms, for which there 
is a single ^He2 ro-vibrational bound state, has been 
studied within the adiabatic hyperspherical representa- 
tion [24, 25]. Parker, et al. have studied the Ne-f-Ne-f-H 
system, calculating the J'^^O^ partial wave contribu- 
tion to recombination and collision induced dissociation 
rates [26]. The studies done by Suno et al. for three he- 
lium atoms [24, 25] and for He-|-He-|- Alkali systems [27] 
have included higher partial wave contributions in order 
to calculate the total recombination rate up to tempera- 
tures of at least 10 mK. These results are relevant to the 
buffer gas cooling technique used in cold and ultracold ex- 
periments [28, 29], where three-body recombination can 
lead to dramatic losses. 

In this paper, we study three-body processes involving 
two hydrogen atoms and one alkali atom. If all atoms are 
spin-stretched, these systems are amenable to calculation 
since the H-l-H interaction has no bound state and X-f-H 
has only a single s-wave bound state for all the alkali 
species, X, considered here [30-33]: '^Li, ^SNa, ^^K, ^^Rb 
and ^•^■^Cs. In thermal Alkali-hydrogen mixtures, the re- 
combination rate for the process H-|-H-|-X^XII+H 
is related to the density of the alkali atoms by 
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where nn and nx are the densities of the H atoms and 
alkali atoms, respectively. Note that there is a 2! reduc- 
tion in the rate if the hydrogen atoms are in a conden- 
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sate [34] . The three-body calculations are thus simplified 
by having only one recombination channel, but the pres- 
ence of sharp avoided crossings in the three-body poten- 
tials makes the calculations a challenge in the adiabatic 
representation. These sharp crossings are between dif- 
ferent families of adiabatic potentials corresponding to 
H+H and X+H. 

We will use atomic units throughout unless specified 
otherwise. It is also convenient to convert energies to 
temperature units by dividing by Boltzmann's constant 
/cb, i.e. 3.17x10-6 a.u.=l K. 



II. METHOD 



This definition of the hypcrangles facilitates the sym- 
metrization of the wave function under exchange of the 
two H atoms. 

After rescaling the three-body wave function ^ as 
?A=i?^/^\I', the three-body Schrodinger equation takes the 
form 
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where V{R, ft) includes all the interactions and is the 
hyperangular momentum operator, defined by taking 
to be the quantization axis for the body-fixed frame [36] 
and expressed as: 



After separating the center-of-mass motion, the rela- 
tive motion of the three particles can be represented by 
the mass-scaled Jacobi vectors and Pi 23 [35] : 
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where ri, r2, and are the lab-frame position vectors 
of the two hydrogen atoms with mass rriH and the alkali 
atom with mass mx, respectively. In the above equa- 
tions, the mass scaling factor d is given by 
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(4) 



The three-body reduced mass is defined as follows to pre- 
serve the phase-space volume element [35]: 



M = 



2mH + mx 



(5) 



In the adiabatic hypersphcrical representation, the hy- 
perradius R, R^=p\2 + P12 3: is the only coordinate with 
the dimension of length and represents the overall size of 
the three-body system. The remaining degrees of free- 
dom, the hyperangles, are represented collectively by fi. 
We use body-frame Delves' coordinates [36] such that 
n= ((/.,0,a,/3,7), with 
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and 6, the angle between the vectors Pi2 and pi 23 such 
that < 9 < TT. The remaining hyperangles are the 
three Euler angles a, /3 and 7 describing the rotation of 
the plane containing the three particles. As a result, the 
interparticle distances rij are determined in terms of the 
internal coordinates {R, 9, (fi) only: 
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The components J = ( Jx, Jy, Jz) are the total orbital an- 
gular momentum operator projected on the body-frame 
axes. 

Since we assume the atoms to be spin-stretched, i.e., in 
the total spin state with the largest magnitude spin pro- 
jection, the relevant Born-Oppenheimer potential surface 
is the lowest quartet surface. We approximate this sur- 
face as a pairwise sum of '^Su two-body potentials: 

y(i?, fi) = uhh(7'12) + vx-ii{r23) + vxiiirsi)- (15) 

The two-body potentials uhh('') and vxair) are shown 
in Fig. 1. At small distances, these potentials are deter- 
mined from ab initio calculations [30, 32, 33] while their 
long-range behavior is determined by the usual dispersion 
potentials [31, 32]. 

Note that for three atoms, a non-additive three-body 
term should be included in V{R,il). This three-body 
term depends on the spatial configuration of the three 
atoms and can be significant for certain configurations. 
The only fully quantum mechanical recombination cal- 
culation for a realistic system that has so far included 
the three-body term found its effect to be negligible [25]. 
But, that work treated He, and He is not very polarizable 
compared to H or the other alkalis. The impact of the 
three-body term for the present systems is thus expected 
to be correspondingly larger. Unfortunately, neither the 
full three-body surface nor the three-body term is avail- 
able for the quartet state of H-l-H-l-X. The importance 
of the three-body term for the systems we are investigat- 
ing, though, can be estimated qualitatively by looking 
at the available three-body terms for the quartet surface 
of identical alkali atoms [37, 38], as they have similar 
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FIG. 1: The potentials for H+H and X+H. Among all 
these combinations, only X+H systems have a single, weakly 
bound molecular state. 



electronic structure. The three-body term is most signif- 
icant when all the atoms arc close together. For Li, the 
three-body potential can make the minimum of the total 
potential about four times deeper than the pair-wise sum 
potential [38]. For heavier alkali atoms, the three-body 
term can change the potential minimum by a factor of 
1.2-1.5 [37]. When the atoms are far apart, a three-body 
dispersion interaction should also be included. The con- 
tribution of this interaction, however, is much smaller. 
For Li, the three-body dispersion interaction is a few 
percent of the pair-wise sum potential, and even smaller 
for heavier alkali atoms [38]. We thus expect that our 
results will change dramatically when the three-body in- 
teractions are included. We know, for instance, that the 
numerical value of the rate can change over a broad range 
for model problems [39]. However, as the first calculation 
for H-l-H-l-X systems, our results can give a sense of the 
order of magnitude for the three-body observables, and 
serve as a starting point for the study of three-body inter- 
actions in such three-body systems. We further note that 
including a three-body term or a full three-body surface 
poses no particular problem for our approach [39]. 

In order to solve Eq. (10), we first expand the three- 
body wave function as 



^1^ = Y,F,e{R)^u{R;^). 



(16) 



u=0 



where fi) are the channel functions obtained as the 

solutions of the adiabatic equation 
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We solve this equation as a function of R, and its eigen- 
values are the adiabatic potentials Uv{R). Therefore, 
upon substitution of -0, Eq. (10) reduces to a set of cou- 



pled ordinary differential equations 
1 (P 



2n dR^ 
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dR 
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F,,e{R)=EF,e{R), (18) 



with non-adiabatic couplings P^^> and Q^^' given by 

d 



dR 



Q^AR)-i'^ 



dR 



(19) 
(20) 



Here, the double brackets denote integration over only 
the hyperangular degrees of freedom. 

In our calculations, the biggest computational burden 
comes from solving the five-dimensional adiabatic equa- 
tion (17). To facilitate its solution, we separate out 
the external degrees of freedom (a, /3, 7) and simulta- 
neously obtain eigenstates of total orbital angular mo- 
mentum and parity by further expanding the adiabatic 
wave functions $jy on the basis of symmetrized Wigner 
D functions [36]: 

J 

^,{R;n) = ^.A'(i?;0,0)^#lf(a,/?,7), (21) 

where H denotes the total parity of the system, and K 
and M denote the projection of the total orbital angular 
momentum on the body-fixed and space-fixed z axes, re- 
spectively. Since molecules XH have only a single s-wave 
state, recombination only happens for the parity-favored 
case, i.e., when H=:(— l)'^. 

Although we are solving for the motion of the nuclei, 
we must require that the total wave function, including 
the electronic degrees of freedom, is antisymmetric under 
exchange of the two protons. For the electronic symme- 
try we are considering, exchanging protons introduces a 
sign change in the electronic wave function. Thus, the 
nuclear part of the wave function must be symmetric un- 
der proton exchange. 

Since we neglect hyperfine interactions, we can cou- 
ple the two protons' spin to give a total spin 1=0,1 
and consider their contribution to recombination inde- 
pendently. For /=1, the spin wave function is symmetric 
under exchange, requiring the spatial wave function to 
also be symmetric. Such spatial symmetry leads to non- 
vanishing K-i at ultracold temperatures [40]. For /=0, 
however, the spin wave function is antisymmetric under 
exchange of the two protons, and the spatial wave func- 
tion is therefore also antisymmetric. For this symmetry, 
K-i vanishes in the zero temperature limit [40]. In the 
present work, we calculate 7^3 for which is the dom- 
inant recombination process for ultracold temperatures. 
For more general cases where the 7=1 state is not pref- 
erentially prepared, our results give only a partial con- 
tribution to Kj, for temperatures beyond the ultracold 
regime. 
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We built the exchange symmetry of the two protons 
into the boundary conditions of the body-frame compo- 
nents Ui,K- Permuting the two protons only affects the 
hyperangles fl: 



Pi26l = n-e. 



(22) 
(23) 



For even parity, the permutation requirements can be 
equivalently expressed as u^k being symmetric about 
0=tt/2 for even K and antisymmetric for odd K. For odd 
parity, u^^k should be antisymmetric for even K and sym- 
metric for odd K. Imposing these boundary conditions, 
we need only solve Eq. (17) in the range < 6 < n/2. 

Asymptotically, i.e., as i? — >■ oo, the adiabatic poten- 
tials with the diagonal couplings included are determined 
by the energies of the break-up components. For the 
atom-molecule channel, the potentials behave like 



(24) 



where the partial angular momentum I is the relative 
orbital angular momentum between the atom and the 
molecule. Since all of the XH systems have only an s- 
wave bound state, l=J. For the three-body break-up 
channels, the potentials behave like 



A(A + 4) 4-15/4 
2/^i?2 ■ 



(25) 



The values of A are non-negative integers determined by 
and the identical particle symmetry [40]. 
Accurate numerical calculations of the three-body ob- 
servables depend largely on the accuracy of the adia- 
batic potentials and channel functions [Eqs. (17) and 
(21)] and, ultimately, on the non- adiabatic couplings 
[Eq. (20)]. By expanding the body- frame components 
Uuk{R]0i4') oh 8- two-dimensional, direct product B- 
spline basis [41], we obtain accurate potentials and cou- 
plings up to i? « 2000 a.u.. Beyond this distance, we 
extrapolate the potentials using the known asymptotic 
expansions [42]. Typically, a [9, 0) mesh of 60x250 gives 
eigenvalues converged to at least eight digits. We have 
found that due to the sharp avoided crossings occurring 
at small i?, a hyperradial grid of about 3000 points is nec- 
essary to accurately resolve most of the abrupt changes 
in the non-adiabatic couplings. Many sharper crossings 
remain, though, that must be traced individually. 



III. THREE-BODY SCATTERING 
OBSERVABLES 



Exa (a.u.) 



a (a.u.) 



Ivdw (a.u. 



H+H 
Li-(-H 
Na+H 
K-FH 
Rb-hH 
Cs-FH 



-1.268x10" 
-3.376x10" 
-7.360x10" 
-2.446x10" 
-1.784x10" 



1.557 
63.71 
43.26 
34.72 
50.24 
56.85 



10.45 
21.50 
22.58 
25.12 
25.92 
27.18 



TABLE I: The two-body bound state energy Exn, scattering 
length a, and van der Waals length l^dw for the H-fH and 
X+H interactions. 



Waals potential —C^lr%j. Therefore, they are character- 
ized by the van der Waals length l^dw = C^f^ijCey^^ [16], 
where fiij is the two-body reduced mass. In Table I, we 
list the bound state energies, the scattering lengths, and 
the van der Waals lengths for all of the two-body poten- 
tials we used. Notice that none the scattering lengths are 
substantially larger than the van der Waals lengths and 
thus the condition for universal behavior (|a| 2> Ivdw) is 
not fulfilled. As a result, we do not expect to observe 
universal physics for these systems. 

After obtaining the potentials and couplings, we solve 
the hyperradial equation (18) using finite elements as de- 
scribed in Ref. [43]. For recombination processes, the to- 
tal recombination rate K3 is the sum over all the partial 
wave contributions [24, 25]: 
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32(2J+ 1) 



I Qjn 



(26) 



where k^y/2jlE and Sj^^ is the scattering matrix ele- 
ment from the initial three-body continuum channel to 
the final atom-molecule channel. From the asymptotic 
form of the three-body entrance channel, the threshold 
behavior of A'3 is determined by the smallest A for that 
symmetry [40] such that: 



cx E'' 



(27) 



In our calculations, we have included the lowest six 
partial-waves J^, implying that A;[^'|j=0, 1, 2, 3, 4 and 5, 
respectively. 

Collision-induced dissociation H+AH^A+H+H is 
the time reversed process of three-body recombination. 
The dissociation rate D3 is defined as [24] 



j,n 



m 

3 



EE 

./,n / 



(2J+l)7r^ 

/^12,3fcl2,3 



I c^n 



(28) 



It is well known [16] that when the scattering length a 
greatly exceeds the characteristic range of the two-body 
interaction, three-body scattering observables are dra- 
matically affected. For the systems we consider here, the 
long-range part of the two-body interaction is the van der 



where ki2,3 = \/2^i2,3{E - Exn) and /.ti2.3 = mH(mH + 
TOjif)/(2mH + mx) is the reduced mass between the H 
atom and the AH molecule. 

Note that the channels that the indices i and / refer 
to are reversed from those in K^. Since the S'-matrix is 



5 



unitary, D3 can be readily calculated once the 5-matrix 
elements for are known. Near the three-body breakup 
threshold where collision-induced dissociation becomes 
energetically possible, D'^^ behaves like [40]: 



a + 2 



(29) 



For atom-molecule collisions, the elastic cross section 
is [25] 



(^2 



j,n 



E 
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(2J-h l)7r 
P 
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cJn 
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(30) 



The threshold behavior of , in contrast to recombina- 
tion, is determined solely by J and follows the standard 
Wigner threshold law. 



\2,/ 



(31) 



The calculation of scattering solutions to Eq. (18) are 
complicated by the sharp avoided crossings in the adi- 
abatic potentials: we typically use 5x10^ hyperradial 
elements distributed as Ri oc from i?=10 a.u. to 
R Ri 2000 a.u.. In the asymptotic region {R > 2000 a.u.), 
the density of elements is fixed to eight elements per 
shortest do Broglie wavelength. To calculate the scat- 
tering observables, we match the numerical solutions to 
the asymptotic analytical solutions at R=1Q^ a.u. for re- 
combination, and at i?=5xl0'^ a.u. for atom- molecule 
collisions. The convergence of the scattering observ- 
ables with respect to the number of adiabatic channels is 
also dramatically affected by the sharp avoided crossings. 
Even the threshold behavior for and requires a 
fairly large number of adiabatic channels for convergence, 
which we take to be from 12 to 25 for all the calculations. 
The resulting K3 and (T2 are converged to at least two dig- 
its for all partial waves, and the three-body bound state 
energies are converged to three digits. 

In our calculations, we have included J^~0^, 1^, 2+, 
3~, 4+ and 5~ for the convergence of the total rates and 
cross sections at high energies. The overall convergence 
of the total total rates and cross sections are converged 
to two digits for E < 200 mK and one digit for 200 mK< 
E <500 mK. 



IV. RESULTS 

A. Three-body recombination rates 

Since the adiabatic hypcrspherical potentials Uv{R) 
arc important in understanding the underlying three- 
body physics involved in the scattering processes, we first 
discuss their behavior. We see that the avoided crossings 
become sharper as we go to heavier alkali atoms. As an 
illustration, in Fig. 2 we show the lowest six adiabatic po- 
tentials U^{R) with jn=0+ for H-hH-hLi and H+H+Cs. 
The potentials for higher partial waves behave similarly 
but become more repulsive as J increases. 
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FIG. 2: (color online) The adiabatic three-body potentials for 
(a) H-f H+Li, and (b) H-|-H4-Cs. The insets show the avoided 
crossings within the energy range we have considered for the 
rates and cross sections. 



To demonstrate the effects of sharp avoided crossings 
on the adiabatic potentials on the non-adiabatic cou- 
plings, we show a key crossing in the insets of Fig. 2. Fig- 
ure 3 shows the corresponding couplings Pi,yi and Qj^jy' 
for H-l-II-l-Cs with J^=Q'^ . It can be seen that when 
a sharp avoided crossing occurs between two potential 
curves, the P^^i and Q^^i coupling those curves show 
sharp spikes at that R. These couplings must be care- 
fully traced out with a dense hyperradial grid in order to 
obtain an accurate solution of Eq. (18). 

In our calculations, the convergence of K-^ depends 
critically on the behavior of the adiabatic potentials. We 
have found that, except for H-l-II-l-Li, all the systems 
have sharp avoided crossings below 1 Kelvin. To cal- 
culate K-i for energies above the crossings like the one 
shown in the inset of Fig. 2(b), both of the potentials 
involved in the crossing need to be included to avoid spu- 
rious resonances. 

The three-body recombination rates for different alkali 
species X are shown in Fig. 4(a)-(f). The correspond- 
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FIG. 3: (color online) The non-adiabatic couplings for 
H+H+Cs with J"=0+. Note that (5i,2 has been multiplied 
by 5000 to have a magnitude similar to Pi, 2- The index v for 
the adiabatic potentials starts from 0. 
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1.3x10"^'' 
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TABLE II: The zero-energy limit of the three-body recombi- 
nation rates K^. 



to more than one digit only for temperatures below 100 
mK. It can be seen that the energy dependence of the 
rates is largely preserved. Further, the thermally aver- 
aged rates for these systems lie close to each other when 
the temperature is beyond 10 mK. In Table II, we list 
the values of K3 for the processes H+H+X^'XH+H in 
the zero-energy limit for reference. 



ing data are available in electronic form [44]. Gener- 
ally, the three-body recombination rates are dominated 
by the = 0+ contribution at ultracold energies, and 
by the j'^ = 1~ and 2+ contributions near the high- 
est energies we have calculated. It can be seen that 
the total three-body recombination rates for H-I-H-I-Na, 
H-|-H+Rb, and H+H+Cs behave similarly and that the 
rates for H+H+Li and H+H+K recombination behave 
differently. In particular, the = 0+ partial rates for 
H+H+K recombination are much smaller than the rates 
for other systems near the zero-energy threshold. The to- 
tal rates for H+H+K recombination are then dominated 
by the = 1~ partial wave contribution for a large en- 
ergy range from about 0.5 mK to 50 mK. Interestingly, 
we have observed that the threshold values of the recom- 
bination rates for different X are ordered by the mag- 
nitude of their non-adiabatic couplings Pqi smd Qoi a-t 
large hyperradii R > 200 a.u. for 0+. This suggests that 
for the present cases 0"*" recombination is dominated by 
inelastic transitions from the lowest continuum channel 
to the atom-moleculc channel at large distances. 

In thermal gases, it is crucial to consider the ther- 
mal distribution of the collisional energies when calcu- 
lating a collision rate. Assuming a Boltzmann distribu- 
tion, we have performed a thermal average of the energy- 
dependent rates. The thermally averaged recombination 
rate {K3) are given by [45]: 

1 r°° 

The results are shown in Fig. 4(f). To perform the 
thermal average, we extrapolate K3 from the lowest en- 
ergy we have calculated to zero energy using the known 
threshold behavior. Since we could not similarly extrap- 
olate to infinite energy for the integral in Eq. (32), the 
integral was limited to the energies we could calculate. 
Consequently, the thermally averaged rates are converged 



B. Collision induced dissociation rates 

Using the simple relation between K3 and -D3 , we have 
also calculated for the same range of E. In Fig. 5, 
we show the thermally averaged collision-induced disso- 
cation rate {D3) as a function of the temperature of the 
XH+H mixture, where {D3) is given by [43] 

The energy E in the integrand is relative to the XH+H 
threshold. That is, it is the XH+H scattering energy, 
and this must be taken into account when evaluating 
Ds{E). It is interesting to note that although dissocia- 
tion is allowed only when the collision energy exceeds the 
molecular binding energy. Fig. 5 shows that in a thermal 
gas dissociation can occur for temperatures well below 
the dissociation threshold. In fact, because we know the 
threshold behavior from Eq. (29) to be D3 a {E — Exh)^ 
for E > ExH (and zero below Exh), we can explicitly 
calculate (D3) below threshold: 

{D3) (X 2Vi(15a;-2)e"^/'= + V^[4+3a:(52:-4)crfc(l/VJ)] 

with X = ksT/ ExH- By contrast, {K3) and ((72) have the 
same threshold behavior as the energy-dependent quan- 
tities (with E replaced by ksT). This formula for {D3) 
is likely valid only for temperatures below i?xH since the 
tail of the thermal distribution starts sampling energies 
outside the threshold regime for higher temperatures, 
making our assumption for the behavior of Dy, invalid. 



C. Atom-molecule elastic cross sections 

As a representative example, we plot the total and par- 
tial cross sections for elastic collisions between H and 
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FIG. 4: (color online) The total three-body recombination rate Kz and the partial rates for H+H+X— ^H+XH, where X 
is (a) Li, (b) Na, (c) K, (d) Rb and (e) Cs. The thermally-averaged total recombination rates are shown in (f). 



KH in Fig. 6. The jn=0+ partial wave contribution 
dominates for collisional energies below 100 mK, beyond 
which the J^=\~ contribution becomes dominant. The 
J^=\~ contribution has a pronounced minimum near 20 
mK, but this feature has only a negligible effect on the to- 
tal cross section. The partial atom-molecule elastic cross 
sections for the other alkali species arc not shown, as their 
energy-dependence is qualitatively the same as shown for 
H-I-KH. Instead, we show in Fig. 7 the thermally averaged 
cross sections for all the alkali species. The thermally av- 



eraged cross sections can be derived from the thermally 
averaged elastic scattering rate, and is given by [43] 



Jo 



(34) 



The total elastic cross sections for all alkali species are 
converged to two digits for all energies. All these data 
are available in electronic form [44]. 

In Table III, we list the values of the elastic cross 
section extrapolated to zero temperature as well as the 
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FIG. 5: (color online) The thermally averaged total collision 
induced dissociation rate (-D3) for H+XH— >-X+H+H. For all 
species, the rates are shown up to 100 mK beyond the three- 
body break-up threshold which is indicated by the vertical 
dashed line. 
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FIG. 6: (color online) The total atom-molecule elastic 
cross section (T2 and the partial cross sections ai^ for 
H-FKH^-H-hKH. 
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FIG. 7: (color online) The thermally averaged total atom- 
molecule elastic cross section (0-2) for H-f XH— >-H-f XH. 





(T2 (cm^) 


47ra^ (cm^) 


au+XH (a.u.) 


H+LiH 


2.2x10"^^ 


1.4x10"^^ 


80 


H+NaH 


8.2x10"^^ 


6.6x10"^^ 


48 


H-hKH 


4.6x10"^^ 


4.2xl0~^^ 


36 


H-fRbH 


1.3x10"^^ 


S.gxlO"^^' 


60 


H+CsH 


1.6x10"^^ 


1.1x10"^^ 


68 



TABLE III: The atom-molecule zero temperature elastic 
cross-section and scattering length between H and XH. 



corresponding values for the atom-molecule scattering 
length an+xK for all the alkali species. Both quanti- 
ties increase with the respective values of the two-body 
scattering lengths (see Table I), or equivalently, the size 
of the molecular state. In fact, we can get an order-of- 
magnitude estimate for the elastic cross-section by sim- 
ply using a2 ~ iira^. These estimated values are also 
shown in Table III. The rough agreement implies that for 
our simple pair-wise sum potential the zero-energy elas- 
tic cross section is mainly determined by the size of the 
XK molecule. By extension, the atom-molecule scatter- 
ing length can be approximated by the X-l-H scattering 
length at the same level of approximation. 



D. Three-body bound state energies 

To complete our study of these systems, we calculate 
the three-body bound states. The three-body energy 
spectra for the H+H+ Alkali systems are very simple due 
to their weakly-interacting nature. In our calculations, 
we have found only one triatomic vibrational bound state 
for ~ 0+ for all the systems. No bound levels 
are found for higher angular momenta. The three-body 
binding energies, relative to the atom-molecule break-up 
threshold, are listed in Table IV. 

To get a sense of the sizes of the triatomic molecules, 
we have also calculated the expectation values of the in- 
teratomic distances (rxn) and (ruH), given by 

poo 

(r-AH)=^/ F,{R)F,,{R){{'S>,\rAii\'^.'))dR, (35) 

where A represents X or H. From these bond lengths, we 
can also calculate the bond angle at the X atom and find 
them to be consistently around 100° for all species. All 
of this geometrical information is included in Table IV. 
From the small binding energies and large bond lengths, 
it is clear that these are very floppy states as is expected 
for van der Waal's molecules. We expect, though, that 
the inclusion of three-body terms in the interaction po- 
tential will tend to bind these states more strongly, reduc- 
ing the bond lengths correspondingly. The three-body 
term may further tend to increase the bond angle towards 
a linear configuration. The three-body term might even 
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Exii2 (a.u.) (rxH> (a.u.) {rnn) (a.u.) Bond angle 



LiH2 9.02x10^* 43 

NaH2 2.58x10"^ 30 

KH2 6.24x10"^ 25 

RbHa 1.95x10"^ 34 

CsH2 1.43x10"^ 38 



65 
46 
37 
52 

59 



98° 
100° 
95° 
100° 

102° 



TABLE IV: The 0^ triatomic bound state energies, expecta- 
tion values of interatomic distances, and bond angles. 



be sufficient to bind additional states, at least for some 
species. 

Finally, we have verified numerically that the triatomic 
binding energies are indeed predominantly determined 
by the lowest adiabatic channel i' ~ 0. Specifically, the 
channel probability J \Fij{R)\^dR for = is beyond 
99% for all the systems. The channel functions Fo{R) 
are shown in Fig. 8. It is interesting to note that, except 
for KH2 , all the triatomic states have a large hyperradial 
extent, reaching values up to a few hundred atomic units 
consistent with the bond lengths listed in Table IV. 



0.25 



02 



0.15 



^ 0.1 



0.05 









KH, 














NaH2 












X^:^^ LiH, 



one alkali atom using a fully quantum mechanical ap- 
proach. Solving the three-body Schrodinger equation in 
the adiabatic hyperspherical representation, we have cal- 
culated the three-body recombination rates and atom- 
molecule elastic cross sections for all the alkali species 
for temperatures up to 0.5 Kelvin. The biggest uncer- 
tainty in our calculations by far is the interaction po- 
tential. Nevertheless, we expect that our results give 
a correct order-of-magnitude estimate of the three-body 
scattering observables. For three-body recombination, 
the lowest three partial waves dominate the total recom- 
bination rates in the energy range we have calculated. 
For the elastic atom-molecule collisions, the cross sec- 
tions are dominated by a single partial wave contribution 
for the energy range in our calculations, which is j'^=0^ 
at lower energies and J^=l~ at higher energies. The 
bound state spectra are very simple, with only one ro- 
vibrational three-body state for each of the alkali species. 

Finally, the difficulty of sharply avoided crossings we 
met at small hyperradius raises an alert for doing adia- 
batic calculations for realistic systems. The complicated 
short-range three-body dynamics can give rise to rapidly 
varying behavior in the adiabatic potentials and the non- 
adiabatic couplings, which makes the adiabatic calcula- 
tions much harder and less reliable. For such cases, a 
diabatic representation of some sort will become neces- 
sary, especially for small distances [46, 47]. 



50 100 150 200 250 300 

R (a.u.) 
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FIG. 8: (color online) The — 0^ hyperradial channel func- 
tion Fo(R) for the triatomic bound states. 



V. SUMMARY 

In this paper, we have studied three-body scattering 
and the bound state spectra for two hydrogen atoms and 
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